FFTS readout for large arrays of Microwave Kinetic Inductance Detectors 
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Microwave Kinetic Inductance Detectors (MKIDs) have great potential for large very sensitive 
detector arrays for use in, for example, sub-mm imaging. Being intrinsically readout in the frequency 
domain, they are particularly suited for frequency domain multiplexing allowing ~1000s of devices to 
be readout with one pair of coaxial cables. However, this moves the complexity of the detector from 
the cryogenics to the warm electronics. We present here the concept and experimental demonstration 
of the use of Fast Fourier Transform Spectrometer (FFTS) readout, showing no deterioration of the 
noise performance compared to low noise analog mixing while allowing high multiplexing ratios. 

PACS numbers: 07.57.Kp 



The Far-Infrared (FIR) and sub-mm wavelength bands 
(10- 1000 /im) contain a plethora of information on the 
evolution and formation of galaxies, stars and planetary 
systems. Future FIR and sub-mm space missions, such as 
SPICA [T] and Millimetron [5] require large arrays of de- 
tectors with an unprecedented high sensitivity, expressed 
in a Noise Equivalent Power (NEP) < 2 x 10~^^W/Hz^/^ 
Also, many ground based telescopes like IRAM [3], 
APEX [4J and future observatories like CCAT [5] require 
very large arrays for faster observation speeds and larger 
instantaneous fields-of-view. 

Microwave Kinetic Inductance Detectors MKIDs |B] 
are a rapidly developing new detector technology 13 |H1 H] 
that has the potential to be used to fabricate very large 
imaging arrays for the FIR, sub-mm [T^, optical and 
x-ray [TT]. The main advantages of MKIDs over other 
detector technologies are: i) ease of fabrication, ii) the in- 
trinsic adaption for frequency domain multiplexing (FD- 
MUX) using microwave readout signal frequencies and 
iii) operation with a wide dynamic range. Currently 
MKIDs have shown sufficient sensitivities in a lab envi- 
ronment for both ground-based and space-based instru- 
ments [12 and there has been an initial demonstration 
at the CSO telescope on Mauna Kea, Hawaii [10,. 

However, FDMUX of MKIDs has yet to be satisfacto- 
rily demonstrated, particularly anywhere near the (power 
and weight) requirements needed for space. Previous 
measurements have been non-scalable using either analog 
mixing to readout one MKID, which has proven useful in 
single detector characterization in a laboratory environ- 
ment, or digital mixing using specialist chips [T3j . Using 
the latter technique a readout of 16 MKIDs has been 
demonstrated using one commercial demultiplexer card. 
However, the technique is hard to scale to very large pixel 
numbers [13]. In this Letter we present a new MKID 
readout scheme based upon frequency division multiplex- 
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ing and a digital Fast Fourier Transform Spectrometer 
(FFTS) OUS]. Digital FFTS systems are presently be- 
ing used as back-end electronics for heterodyne mixers on 
several ground based telescopes [11] . They can process a 
real time data stream with a bandwidth of up to 1.8 GHz 
into a frequency spectrum of ~ 8000 points. We propose 
here to use the FFTS to read out a set of single frequency 
probe signals that have passed through a chip contain- 
ing one MKID detector per probe signal. The proposed 
readout is able to read out ~ 1000 pixels using 1 probe 
signal generator and 1 FFTS board. To demonstrate 
the principle, we describe an experiment where we have 
read out 8 MKIDs simultaneously and we have measured 
both the dark detector NEP and the response to an op- 
tical signal for 8 MKIDS. These results are compared to 
a conventional analog single MKID readout scheme and 
show that the FFTS based multiplexed readout does not 
deteriorate the intrinsic system noise of our setup, which 
is limited by the first stage cryogenic amplifier. 

MKIDs are superconducting pair breaking detectors 
that sense the change in the complex surface impedance 
of a thin superconducting film due to radiation ab- 
sorption with a (sky) frequency ¥rad > 2A/h, which 
is ~ 80 GHz for aluminum. An MKID consists of a 
thin superconducting film that is incorporated in a res- 
onance circuit which is either capacitively [6] or induc- 
tively [TT] coupled to a through line. Changes in the sur- 
face impedance of the film are converted to changes in 
resonator quality factor and resonance frequency. These 
changes can be read out by measuring the phase and am- 
plitude modulation of a probe tone at a probe frequency 
equal or close to the resonator resonance frequency Fq, 
which is typically a few GHz. 

MKIDs take advantage of the fact that a supercon- 
ductor at Tc/10 has negligible losses at the probe fre- 
quency /q. This enables very low 3 dB bandwidth Sf 
resonators with high Q factors fo/6f ^10^. Hence, close 
packing of the resonators in (readout) frequency space is 
only limited by manufacturing tolerances (~1 MHz at 
/o ^3 GHz). For example, using conventional one oc- 
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FIG. 1: The circuit diagram for test in this experiment. The 
readout signal is generated in software from the known MKID 
resonance frequencies in the IF (0.1 -300 MHz) and converted 
to analog by the DAC of the AWG card; the signal is then 
mixed with signal generator 1 (SGI) up to ~4-4.3 GHz using 
an upconverter (inset lower right), passes through the MKIDs, 
amplified and down converted back to the IF band. A last 
amplification step follows before digitization and processing 
by the FFTS. For comparison, by using a microwave switch 
we can switch from the FFTS readout to a conventional single 
pixel readout scheme using an IQ mixer and a separate LO 
(SG2) at the resonance frequency of one of the MKIDs. 



tave bandwidth amplifier (4-8 GHz) one can read out 
~ 4000 MKIDs using just one pair of coaxial cables. 

To demonstrate the FFTS based readout of MKIDs 
we use the readout scheme shown in Fig. [l] to readout 8 
MKIDs. A set of 8 Intermediate Frequency (IF) probe 
tones in a band from 0.1 -300 MHz are generated by a 
1 GSPS (1 X 10^ samples per second) commercial arbi- 
trary waveform generator (AWG) with a 12-bit Digital- 
to- Analog Converter (DAC). Subsequentially they are 
mixed to the required readout at Radio Frequency (RF) 
of 4 -4.3 GHz by a commercial single sideband upcon- 
verter using an Local Oscillator (LO) frequency of 4 GHz 
from SGI (example tones shown in the inset Fig.[lJ. The 
resulting RF probe signal consists of a set of RF fre- 
quency probe tones, where each tone corresponds to one 
MKID while additional blind tones can be added to mea- 
sure the system noise. The LO of the mixer is a low noise 
synthesizer which also has an 1 GHz clock output used to 
clock the DAC and FFTS. The tones are run at as high 
a level as possible to use the full scale of the DAC and 
minimize the level of the LO leakage into the RF signal. 
Subsequently, the signals are attenuated using microwave 
attenuators at room temperature, 4 K and 1 K to the re- 
quired level MKID readout power of ~ —90 dBm per 
MKID. The signals pass through the MKID chip that is 
cooled to a temperature of 100 mK. The tones are mod- 
ulated by the MKIDs, amplified at 4 K by a low noise 
amplifier (T„ = 4 K) and passed out of the cryostat. Here 
a microwave switch enables the use of the FFTS readout 
or a conventional single tone readout (IQ readout) [12 of 



one of the 8 multiplexed MKIDs. For the single tone IQ 
readout, the RF probe signal is passed to an IQ (quadra- 
ture) mixer where the LO is generated by a second signal 
generator (SG2) tuned to the exact probe tone frequency 
of one of the resonators. The two signal generators, DAC 
card and FFTS board are phase locked using a 10 MHz 
reference clock to reduce the relative clock jitter. By 
changing the SG2 frequency we can switch between all 8 
probe tones. The IQ mixer outputs the In-phase (I) and 
Quadrature-phase (Q) signals of the tone of the selected 
MKID at baseband that is then digitized using a dual 
channel 16 bit ADC at 2x10^ samples/s. 

For the multiplexed FFTS readout, the RF probe sig- 
nals go to a downconverter which is driven by the same 
LO as the upconverter (SGI). As a result the RF probe 
tones are converted back to the exact same IF frequen- 
cies as the original IF tones. A final amplification is used 
to compensate for conversion loss at the mixer to effi- 
ciently use the full scale deffection of the FFTS Analog 
to Digital Converter (ADC), which digitizes the signal 
at 8 bit and 2 GSPS. The FFTS performs a continu- 
ous N=8192 bin (point) complex Fast Fourier Transform 
(FFT) outputting a N bin power spectrum for every block 
of N points of time domain data. The outputted 8192 bin 
spectrum has a binwidth of 1 GHz/8192=:122 kHz, where 
the value of each bin is the summed spectral power in the 
bin from that block. By comparing the spectrum from 
different blocks of time domain data we can follow the 
time evolution of the spectral power in each bin. By ar- 
ranging that the probe frequencies occupy different bins 
this one board will be equivalent of 8192 digital down- 
converters [T5 . A selection of the bins have a probe tone 
in them allowing the time evolution of the probe tones to 
be simultaneously followed at the FFT rate of 122 kHz, 
intrinsically multiplexing the readout of the tones in the 
frequency domain. Note that we measure now the mod- 
ulation of the probe tone power which is a measure of 
the power transmission of the MKIDs, |S'2iP- The signal 
bandwidth is then reduced by averaging multiple spectra 
for 20 ms or more. This final integration bandwidth sets 
the noise bandwidth of the readout, not the bin band- 
width. 

The devices used in the experimented are coplanar 
waveguide quarter wavelength resonators [6l [12] capaci- 
tively coupled to the through line. The devices are made 
from a single 100 nm thick Al film thermally evaporated 
on a high purity silicon substrate. The measured chip 
has 25 usable MKIDs at slightly different resonance fre- 
quencies. 8 MKIDs with Q ~ 150000 were chosen for the 
experiment. 

In the first experiment we measure the response of 8 
probed MKIDs to a red LED illumination of the chip with 
optical photons, increasing the LED current in discrete 
steps [12]. The insets of Fig. [2] show the measured fre- 
quency dependent transmission and probe frequencies for 
one MKID. Fig. ^a) shows the multiplexed (i.e. FFTS 
readout) response of 8 MKIDs with the probe tones set 
on resonance when there was no LED illumination. We 
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FIG. 2: (a): Measured response of SMKIDs in time to a op- 
tical step ramp from a red LED. The initial probe frequency 
was on resonance, the LED was stepped linearly in current, 
(a inset): transmission curve versus measurement fre- 

quency measured by the FFTS technique of a MKID around 
the resonance frequency. The frequency was swept by chang- 
ing the LO frequency of the mixers, keeping the DAC output 
constant. Indicated are the approximate probe point and the 
off resonance frequencies, (b): Power spectral density of noise 
compared to carrier in the signal amplitude at off resonance 
(lower) and the probe point which was 8 kHz from resonance 
(upper), measured using the FFTS readout and compared to 
measurements using an IQ mixer. 



clearly observe the expected step wise increase in 15211^. 
Non uniform illumination of the chip is the reason of ob- 
served response difference. 

To fully demonstrate the readout we need to study 
not only the MKID response but also we need to study 
the noise properties of the new readout technology. To 
do this, we sample |S'2iP of each probe signal using the 
FFTS at a rate of one measurement every 20 ms for 10 s. 
In this experiment we have used probe tone frequencies 
at ^ 0.6(5/ away from the resonance frequency. The mea- 
surement is repeated by tuning the probe signal 1 MHz 
off resonance, effectively measuring the system noise not 
influenced by the MKID. This data is used to calculate 
the power spectral density with respect to the carrier am- 
plitude, the result is shown for one representative of the 
8 MKIDs by the bold lines in the figure, both for 8 kHz 
and 1 MHz off resonance. Subsequently, we measured the 
noise spectrum for each MKID separately in time both 
at 8 kHz and 1 MHz off resonance by using our conven- 
tional IQ setup. The result is shown by the thin lines in 
the Fig. (2) 

We observe several features. First, we observe that 
there is no measurable difference between the FFTS read- 



out traces and the IQ readout traces. Since the latter is 
limited by the setup, we conclude that the FFTS readout 
system does not limit the system noise. Secondly, both 
setups show a large excess noise when probing the MKID 
8 kHz off resonance, which is roUed-off at ^15 kHz. This 
is the magnitude representation of the MKIDs excess fre- 
quency noise [7] , rolled off at the ring response frequency 
of the resonator (bandwidth Sf = fo/{2Q) ~ 15 kHz for 
the MKIDs presented here) . Above this roll-off the noise 
returns to the system noise floor. Thirdly, we observe 
for both datasets a much lower, almost white noise spec- 
trum when probing the MKIDs at 1 MHz off resonance. 
This noise level is consistent with the system noise tem- 
perature. The small difference at the highest frequencies 
between the 1 MHz and 8 kHz off resonance data is identi- 
cal to the difference in ^21 between the 8 kHz and 1 MHz 
off resonance transmission, a direct consequence off our 
normalization to the carrier amplitude. 

From the measured noise trace at 8 kHz, tempera- 
ture dependent response of the MKID and in combina- 
tion with a separate measurement of the quasiparticle 
lifetime (~ 1 ms) we can estimate a dark NEP [20j of 
- 1 X lO-i^W/Hzi/2 at 20 Hz. 

The propose readout of 8 pixels can easily be expanded 
to much more pixels. With the existing boards the fun- 
damental limit is given by the 8192 bins of the FFTS. 
However, for so many pixels careful control is required of 
the peak voltage at each stage in the readout chain [TS] 
to prevent clipping at the ADC and to control intermod- 
ulation products from the non-linearity of the amplifiers 
when overdriven. For a signal consisting of n tones with 
identical rms power P, the peak power can be approxi- 
mated by Ppeak PnC , where n is the number of tones 
and C the Crest factor, where C — Ppcak/nP ~ 25 using 
random tone phase. The consequence is that the power 
gain available between the MKIDs and ADC cx n. Since 
the contribution of the ADC to the system noise (referred 
to the input of the MKIDs) is dependent on the system 
gain, increasing the amount of probe signals requires a 
lower noise from the ADC |16j which can be expressed in 
terms of the bit noise. Effective Number Of Bits (ENOB). 
For example, consider a 1000 MKID pixel array with a 
Q factor of 50000 and probe signal power per MKID of 
-85 dBm. This would require an ADC with - 9 ENOB 
to give the effective noise at the chip level identical to 
the low noise amplifier noise floor. In the case of phase 
or 15*211^ readout, the system noise is higher due to the 
excess device noise, so relaxing the requirements for the 
ADC in the example given to about 7 ENOB. 

The readout described in this letter gives out |S'2iP 
response of a MKID. A readout in the complex plane 
is preferred [121 giving the phase and amplitude change 
of the probe signals, enabling either a more linear re- 
sponse (phase readout) or a higher sensitivity for low 
background applications using high Q MKIDs from both 
the phase and amplitude of 5*21. The current FFTS does 
internally have the complex information from which the 
complex transmission can in the future be extracted. 
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In conclusion, we present a MKID readout concept 
which will facilitate the application of MKIDs to a wide 
range of applications. The electronics, as presented here, 
can already handle 100 pixel readouts for ground-based 
mm/sub-mm astronomy, easily expandable with wider 
bandwidth tone generation schemes or closer readout 
packing of the MKIDs. This can then be expanded with 



a phase/amplitude readout to reach the more stringent 
requirements for very low NEP applications like space 
born sub-mm spectroscopy, for example SPICA [T]. 
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